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Abstract: Within the year 2007, MAGIC will be upgraded to a two telescope system at La Palma. Its 
main goal is to improve the sensitivity in the stereoscopic/coincident operational mode. At the same time 
it will lower the analysis threshold of the currently running single MAGIC telescope. Results from the 
Monte Carlo simulations of this system will be discussed. A comparison of the two telescope system 
with the performance of one single telescope will be shown in terms of sensitivity, angular resolution and 
energy resolution. 



Introduction 

The MAGIC telescope, currently the largest (17 m 
diameter mirror) single dish imaging atmospheric 
Cherenkov telescope, has been in scientific oper- 
ation since summer 2004. A second 17 m tele- 
scope equipped with advanced photodetectors and 
ultra-fast readout is under construction and is ex- 
pected to be ready in 2007 [1]. MAGIC-II, the two 
telescope system, is designed to lower the energy 
threshold and simultaneously achieve a higher sen- 
sitivity in stereoscopic mode. 

In order to determine the optimal baseline distance 
between the first and the second telescope and 
estimate the performance of the system, detailed 
Monte Carlo simulations have been carried out [2]. 
The Monte Carlo simulation of the MAGIC-II sys- 
tem is divided into three stages. The CORSIKA [3] 
program simulates the air showers initiated by ei- 
ther high energy gammas or hadrons. In this simu- 
lation we have used the CORSIKA version 6.019, 
the EGS4 code for electromagnetic shower gener- 
ation and VENUS and GHEISHA for high and low 
energy hadronic interactions respectively. New at- 



mospheric models have been introduced on the ba- 
sis of studies of total mass density as a function of 
the height. The second stage of the simulation. Re- 
flector program, accounts for the Cherenkov light 
absorption and scattering in the atmosphere and 
then performs the reflection of the surviving pho- 
tons on the mirror dish. Finally, the Camera pro- 
gram simulates the behaviour of the MAGIC pho- 
tomultipliers, trigger system and data acquisition 
electronics. Pulse shapes, noise levels and gain 
fluctuations obtained from the MAGIC data have 
been implemented in the simulation. 

For the present study a total of 1.14 x 10^ pro- 
tons between 30 GeV and 30 TeV have been pro- 
duced, as well as and 2.0 x 10^ gammas between 
10 GeV and 20 TeV. The energy distribution of pri- 
mary gamma rays is a pure power law with a spec- 
tral index of -2.6, whereas charged primaries fol- 
low the power law of -2.78. The telescope point- 
ing direction is 20° in zenith, with the directions 
of protons scattered isotropically within a 5° semi- 
aperture cone around the telescope axis. Maximum 
impact parameters of 350 and 450 m have been 
simulated for gammas and hadrons respectively. 



MAGIC-II Monte Carlo simulation 



Analysis of stereo events 

The two telescopes in the MAGIC-II system can 
be independently operated by observing two dif- 
ferent sources or sky regions. However, the best 
performance of the system is achieved with the si- 
multaneous observation of air showers by the two 
telescopes. The stereoscopic observation mode al- 
lows a more precise reconstruction of the shower 
parameters as well as a stronger suppression of the 
hadronic showers and other background events. 

The analysis of stereoscopic events is performed 
by individually analyzing the images from the two 
telescopes. A set of parameters (Hillas parame- 
ters [4]) is obtained from each image and they are 
combined to obtain the shower parameters. Only 
showers triggering both telescopes are considered 
under the stereo analysis. The images are com- 
bined following the first algorithm in [5]. The in- 
tersection point of the two major axis of the ellipses 
recorded in the telescope cameras, provides the lo- 
cation of the source of a particular shower (fig- 
ure 1). The 6'^ parameter is defined as the square 
of the angular distance from the real source image 
in the camera and the reconstructed one for each 
event. The location of the shower core on ground 
is obtained by intersecting the image axes from the 
telescope positions on the ground. In addition, the 
height of the shower maximum (Hmax) can also be 
obtained. Having only two telescopes, the quality 
of the reconstruction of these parameters depends 
on the amplitude (Size) of the recorded images 
and the angle between the axes. 

Image analysis is based on the analysis of the sec- 
ond moments of the images recorded by the tele- 
scopes cameras. Among the parameters obtained. 
Width and Length (second moments of the light 
distribution along the major and minor axes of the 
image) are the more useful ones for signal (gam- 
mas) / background (hadrons) separation. Width and 
Length depend on the distance from the shower 
core to the telescope, hence, a correct estimation of 
the impact parameter is required to properly eval- 
uate these parameters. With a single telescope, the 
observer can not easily resolve the ambiguity be- 
tween a close by, low energy shower and a dis- 
tant, high energy one. With a second telescope, 
in most cases the ambiguity disappears because of 
the stereoscopic vision of the showers. 




Figure 1: Stereoscopic principle. The incoming 
direction of the shower and core distance are ob- 
tained by combining the information of the images 
on both telescopes. 

In order to combine the parameters from both im- 
ages, we compute the Mean Scaled Width(MSW) 
and Length(MSL) parameters. These new param- 
eters are obtained by subtracting the mean and 
dividing by the RMS of the parameter distribu- 
tion (as a function of size) for Monte Carlo gam- 
mas. The new distributions have a mean value 
of and RMS of 1 for gamma showers and are 
broader for proton showers. The distribution of 
these parameters for gamma and proton showers 
are shown in figure 2. A comparison between 
the single telescope case (MAGIC-I) and two tele- 
scopes case (MAGIC-II) is also shown. A higher 
gamma/hadron separation is achieved when com- 
bining the information from the two telescopes, 
compared with the single telescope case. 

In our analysis, the Random Forest [6] (RF) tech- 
nique is used for gamma/hadron separation. RF is 
a multidimensional classification tool that, in this 
case, it is used to determine an average probability 
of an event to be a hadron induced shower. A set 
of MC gamma ray events and MC proton events is 
used to train the RF. After training, the test sam- 
ples can be classified by RF, providing a param- 
eter (called Hadronness) distributed between 
and 1. Values closer to mean that the event is 
more gamma-like and values closer to 1 mean that 
the event is more hadron-like. 

For this study, the parameters that have been used 
in Random Forest are: average amplitude {Size), 
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Figure 2: Comparison of the scaled Width (left) 
and Length (right) parameters between MAGIC-I 
and MAGIC-II. 

average number of islands (Nisiand) \ core dis- 
tance, Hmax, MSW and MSL. This set of pa- 
rameters, however, should not be considered as op- 
timal and some improvement could be expected 
with an optimized parameter selection. 

In the analysis, gamma/hadron separation is based 
on the Hadronness and 0'^ parameter is used to 
extract the signal events. 

Finally, energy reconstruction is based on lookup- 
tables where Monte Carlo energy of gamma show- 
ers is tabulated as a function of shower impact pa- 
rameter, Hmax and Size. For each telescope, a 
reconstructed energy is obtained by interpolation. 
The shower energy is obtained as the average of 
the values from both telescopes. 

Simulation results 

Monte Carlo simulations show that the sensitivity 
does not vary dramatically with the distance be- 
tween telescopes, the optimal value being around 
90 m [2]. The sensitivity is defined as "integral 
flux resulting in gamma excess events, in 50 hours 
of observation, equals to 5 times the standard devi- 
ation of the background". The effective area at trig- 
ger level of the MAGIC-II system is smaller than 
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Figure 3: Comparison of the angular resolution be- 
tween MAGIC-I (1 telescope) and MAGIC-II (two 
telescopes). Sigma of the fit to a 2-dimensional 
gaussian is also shown for MAGIC-II. 

that of any of the two telescopes because of the co- 
incidence requirement. However, at energies above 
300 GeV matches the effective area of MAGIC-I 
which is slightly lower than that of the new sec- 
ond telescope because this is equipped with higher 
quantum efficiency photomultipliers. The coinci- 
dence requirement provides a higher background 
rejection. As a result a better sensitivity below 100 
GeV and also a reduction of the analysis threshold 
is achieved. 

For a single telescope, the angular resolution is es- 
timated using a modified parametrisation of the so 
called DISP method [7]. The discrimination of 
the shower head and tail relies on the shape of the 
image (asymmetry along the major axis). This of- 
ten results in a wrong head-tail assignment that de- 
grades the angular resolution. With two telescopes, 
this drawback is easily overcome since the source 
direction is obtained as the intersection of major 
axes of the images in the camera. The angular reso- 
lution, here defined as the angle within which 50% 
of the reconstructed gammas from a point source 
would be contained, as a function of gamma ray 
energy is shown in figure 3. The improvement in 
angular resolution from 1 to 2 telescopes is clearly 
seen. 

The stereoscopic analysis also results in a better 
energy reconstruction due to better reconstruction 



1 . an island is defined as any cluster of 2 or more pix- 
els surviving the image cleaning. 
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Figure 4: Energy resolution of the MAGIC-II (two 
telescopes) system compared with MAGIC-I. 
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Figure 5: MAGIC-II system sensitivity compared 
with the Crab flux and other existing experiments 
(MAGIC-I, HESS and VERITAS). 

of the shower axis and also a double sampling of 
the light pool. The energy resolution for gammas 
as a function of primary energy is shown in fig- 
ure 4. For comparison, the energy resolution of 
MAGIC-I is also shown. An energy resolution 
for gammas better than 20% is achieved above 50 
GeV. 

The sensitivity estimate for MAGIC-II is shown in 
figure 5. The flux sensitivity of the 2-telescope sys- 
tem is between 2 and 3 times better than that of a 
single telescope (MAGIC-I) and it is significantly 
improved below 100 GeV. The MAGIC-II system 
can achieve a sensitivity of 1% Crab in 50 hours 
above 150 GeV. 



Conclusions 

The MAGIC-II system of 2 telescopes will per- 
form observations in stereoscopic mode, allowing 
a more precise reconstruction of the showers and 
a significant reduction of backgrounds below 100 
GeV. This will make possible an improved angu- 
lar and energy resolution as well as a reduction of 
the analysis threshold. All together it will increase 
the current sensitivity of the instrument by a factor 
between 2 and 3 at different energies. 
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